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THE  LATERAL  BUCKLING  OF  A  STRAIGHT  OR 
CURVED  BEAM  SUBJECTED  TO  PURE  BENDING* 

by 

F.  E.  Vanslager+ 

Department  of  the  Aerospace  and  Mechanical  Engineering  Sciences 
University  of  California,  San  Diego 
La  Jolla,  California 


ABSTRACT 

This  paper  presents  a  new  analysis  of  the  classical  problem  of 
the  lateral  buckling  under  pure  bending  of  a  curved  beam  with  circular 
axis.  Although  a  solution  to  this  problem,  and  to  the  corresponding 
straight  beam  problem,  is  well-accepted  in  the  literature,  the  present 
analysis  gives  buckling, moments  vrhich  are  significantly  different  than 
the^usually  accepted  results.  The  difference  arises  from  an  application 
of  new  boundary  conditions  tnat  are  shown  to  be  the  only  boundary  con¬ 
ditions  consistent  with  the  equilibrium  of  the  buckled  beam  under  a 
pure  bending  moment.  ThiV  new  analysis  provides  generally  higher 
bending  moments.  A  short  experimental  program  was  performed,  and 
the  results  tend  to  confirm  the  new  analysis. 


The  results  presented  in  this  paper  were  obtained  in  the  course  of 
research  sponsored  under  Contract  No.  N00014-67-A-0109-0003 
Task  NR  064-496  by  the  Office  of  Naval  Research,  Washington,  D.C. 
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arbitrary  constants 
arbitrary  constants 

torsional  rigidity  a!,  out  the  centroidal  axis 
Young's  modulus 

second  moment  of  the  area  about  the  x  and  y  axes 

an  eigenvalue  of  eq.  (3)  or~(23) 

an  eigenvalue  of  eq.  (3) 

length  o*  a  straight  beam 

bending  moment  vector 

components  of  the  bending  moment  vector 
an  integer 
axial  force 
mean  radius 

coordinate  along  the  centroidal  axis 
displacement  in  the  x  direction 
work  of  deformation 
displacement  in  the  y  direction 
potential  energy 
work  of  external  forces 
coordinate  in  the  radial  direction 
coordinate  transverse  to  the  beam 
beam  angle 
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1.  INTRODUCTION 


This  paper  is  concerned  with  the  well-known  problem  of  the 
lateral  buckling  under  pure  bending  of  a  straight  beam  or  of  a  curved 
beam  with  a  circular  axis.  A  sketch  of  a  curved  beam  of  rectangular 
cross  section  with  mean  radius  R  and  opening  angle  a  is  shown  in 
Figure  1.  It  is  assumed  that  the  beam  is  loaded  by  moments  M  acting 


at  s  =  0  and  s  =  aR  and  is  otherwise  unloaded.  In  order  to  have  a  con¬ 
servative  system,  it  is  also  assumed  that  the  direction  of  the  applied 
moment  vector  does  not  change  during  the  buckling  process.  When  the 
applied  moment  M  reaches  a  critical  value,  the  planar  configuration 
becomes  unstable,  and  the  beam  begins  to  deform  laterally.  Section 
A-A  shows  the  displacement  of  a  typical  cross  section  of  the  buckled 
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beam;  the  initial  position  is  indicated  by  the  dashed  rectangle. 

The  classical  solution  to  this  problem,  obtained  by  S.  Timoshenko 
and  published  in  Bull.  Polytech.  Inst.  Kiev  in  19’ 0,.  is  widely  accepted, 
viz  [l],  [2]*,  and,  tc  the  author's  knowledge,  it  has  never  beendisputed. 
In  the  first  part  of  the  present  paper,  a  different  solution  is  obtained  for 
both  the  circularly  curved  beam  and  the  straight  beam.  The  last  part  of 
the  paper  describes  an  experimental  verification  of  the  new  solutions  as 
opposed  to  the  classical  solutions.  The  discrepancy  between  the  two 
solutions  is  shown  to  be  due  to  the  boundary  conditions  prescribed  on  the 


loaded  ends  s  =  0  and  s  =  aR.  The  only  conditions  that  are  consistent 
withuthe  basic  equations  [l]  are  shown  to  be  those  of  zero  geometric 


constraint  —  i.  e. ,  "free  edges",  and  the  solution  for  these  conditions  is 
derived  in  section  3  below. 


2.  BASIC  EQUATIONS 

The  linearized  strain-displacement  relations  for  the  transverse 
curvature  x^,  the  in-plane  curvature  change  x^.,  and  the  unit  twist  |3 
are  [1],  [3], 


eg  d3  v  _u_  d3u 

x  R  ds2  ’  *y  R2  ds2 


where  u  and  v  are  linear  displacement  components  and  is  a 
rotation  as  shown  in  Figure  1.  The  stress -strain  relations  are  given  by 


X  =  ~ 
x  El 


xy  =  Ei 


Numbers  in  brackets  refer  to  references  in  the  List  of  References. 
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3 


(6) 


< 1 5,  =  a,  sinks  where  k  =  — —  and  n  =  0,  ±1,  ±2,  ... 
k  k  aR 

or,  alternatively,  to  apply  the  "built- in"  boundary  conditions  that 
d<Pk 

**= — -  =  0  at  s  =  0  and  s  =  ali,  which  leads  to 
da 

=  b^  cos  ks  where  k  »  ~~  and  n  =  0,  ±1,  ±2,  ...  .  (7) 

However,  an  inspection  of  equation  (4)  shows  that  both  of  these 
solutions  lead  to  a  finite  positive  buckling  moment  for  an  infinitely 
long  beam,  since  equations  (4),  (6),  and,  (7)  show  that  k  -*  0 

o 

(M  -*  C/R  or  EI^/R)  as  a ®> .  Another  objection  is  that,  as  shown 
by.  equations  (3)  and  (4),  the  minimum  buckling  moment  M  =  0  occurs 
when  k  =  —  ,  and  corresponds  to  a  rigid  body  rotation  about  a  diameter*. 
Therefore,  since  equations  (6)  and  (7)  would  give  k  =  4*  for  a  =  it  and 
n  =  1,  the  buckling,  moment  would  change  from  a  small  negative  value 
to  a  small  positive  value  as  the  beam  angle  a  passes  through  tt.  A 
plot  of  the  neg?/rive  buckling  moment  vs.  beam  angle  a  would,  there¬ 
fore,  show  a  discontinuous  jump  to  the  next  higher  mode  at  a  =  ir.  And, 

finally,  since  the  rigid  body  rotation  solution  with  k  =  —  should  be  in- 

R 

dependent  of  the  angle  a,  and  since  equations  (6)  and  (7)  show  that 
k  =  ~  occurs  only  for  discrete  values* of  a.,  we  conclude  that  equations 
(6)  and  (7)  are  probably  not  the  correct  solutions  to  the  problem. 

What,  then,  is  wrong  with  this  classical  approach?  For  an 


*  One  could  equally  well  take  k  =  -—  as  the  fundamental  rigid  body 
rotation  but  this  has  some  notational  disadvantages. 
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answer  we  must  take  a  closer  look  at  the  statement  of  the  problem.  It 
is  simply  the  question  of  determining  the  magnitude  of  the  equal  and 
oppositely  directed  end  moments  M,  at  which  the  team  begins  to  later¬ 
ally  deform  in  some  unknown  fas  hi  on.  Equation  (3)  is,  therefore, 
derived  by  assuming  that  the  deformed  beam  is  in  equilibrium  with  an 
end  moment  M  which  moves  with  the  beam,  but  whose  direction 
remains  fixed  in  space.  Since  no  geometric  boundary  conditions  are 
implied  in  the  derivation  of  equation  (3),  the  ends  are  essentially  "free", 
and  any  solution  of  equation  (3)  is  potentially  a  solution  of  the  problem. 

In  general,  however,  arbitrary  boundary  conditions  will  give  rise  to 
shearing  forces  on  the  boundaries  which  were  not  included  in  the  deriv¬ 
ation  of  equation  (3).  Therefore,  additional  principles  must  be  invoked 
in  order  to  determine  what  are  the  correct  boundary  conditions  for  the 
problem.  It  turns  out  that,  since  k  =  ^  is  a  solution  to  the  problem, 
it  is  sufficient  to  require  that  the  principle  of  stationary  potential 
energy  be  satisfied;  i.  e. ,  in  the  vicinity  of  the  true  solution,  t‘*»e  vfoi'k 
of  deformation  must  be  equal  to  the  work  of  the  applied  moment  M. 

In  evaluating  the  potential  energy  for  the  beam,  we  will  assume 
that  1^  is  very  large  in  comparison  to  and  C.  The  work  of  defor¬ 
mation  U  will  then  be  given  by 


U 


El  aR  r  paR 

~  f  x  3  ds  +  ~  f  p3ds 
2  J0  x  2  JQ 


or,  by  using  equation?  (1)  and  (2), 
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there  is  a  small  change  in  curvature  of  amount 

1  I  _  I  __1 _ _A  R(cos  0-1)  _  Q 

R*p“R“R-A~“R2=  R2  ~  "  2R 


Integrating  these  curvature  changes  along  the  length  of  the  beam  to  get 
the  total  angle  change  and  multiplying  by  the  applied  moment  M,  we 
get  that  the  work  W2  due  to  the  "rotation"  of  the  ends  is  given  by 


W_  =  ^  .  M 


\EI 


2R 


(10) 


Subtracting  W_  and  W  from  the  work  of  deformation  U,  we  find, 

1  2 

after  some  rearrangement,  that  the  potential  energy  V  can  be  written 
as 


V  = 


M 

2EI 

x 


aR  ,  3 

J#  (f) ds-  (n> 


Using  the  definition  of  k  given  by  equation  (4),  this  can  also  be 
written  as 


In  this  form  it  can  be  readily  seen  that  the  Euler -Lagrange  equation 
corresponding  to  the  minimization  of  the  potential  energy  V  is  the 
equilibrium  equation  (3).  This  is  as  it  should  be,  indicating  that  the 
potential  energy  has  been  correctly  derived. 

We  can  now  apply  the  principle  of  stationary  potential  energy. 
We  take  a  varied  state  in  the  form 
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(13) 


<p*  =  <Pk,  +  € 

where  is  the  assumed  "true"  solution  to  equation  (3)  under  theapplied 
moment  determined  by  k  from  equation  (4),  is  the  "true"  solution 

V 

of  equation  (3)  with  k  replaced  by  l,  and  e  is  an  arbitrarily  small 
constant.  Substituting  (13)  into  ( 1 2)  and  requiring  that  the  potential 
energy  be  stationary  at  the  true  solution  (p  ,  we  obtain 

K 


Kf 

MC 


Cn 

)  aviol 

d€ 


=  k2 


OR  OCR  d<p d<p, 

J0  V-id3-J0  (-S-)(dT)d8  =  °' 


(14) 


We  now  take  another  varied  state  in  the  form 


«**=  H  +  ^k 


(15) 


where  now  the  applied  moment  is  determined  by  l,  and  the  roles  of 
ip  and  <p  have  been  reversed.  Then,  requiring  the  potential  energy 

JK 

to  be  stationary  at  the  true  solution  <p  ,  we  obtain 


Thus, 


if,  k2  £  -t2  ,  we  obtain  from  equations  (14)  and  (16) 


(17) 


These  orthogonality  relations  between  the  eigenvectors  are  analogous 
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to  the  orthogonality  relations  for  the  buckling  of  a  straight  beam  under 
an  axial  force  (see  the  Appendix).  By  integrating  equations  (17)  by 
parts  and  using  equation (3),  one  can  obtain  the  natural  homogeneous 


boundary  conditions  on  the  These  are 


d<^_,aR 


=  0 


and 


(18) 


r  “’’’in0*  „ 

[^-arj0  =  0  • 


The  orthogonality  conditions  are  sufficient  to  determine  the 
eigenvalues  k  uniquely  when  it  is  remembered  that  the  rigid  body 
rotation  solution  must  be  an  allowed  solution  to  the  problem.  In  order 
to  emphasize  that  k  =  ~  is  a  solution  of  this  buckling  problem,  we 
will  rewrite  equations  (5)  in  terms  of  odd  and  even  parts,  with  this 
rigid  body  rotation  solution  explicitly  included  as  follows: 

=  AHr  -  2  )  +  V03  (I  '  I  )  '  (19) 

«>k  =  Aksink(s  -  T)  +  Bkcosk(s  -v)  ' 


When  the  orthogonality  conditions  (17)  are  applied  to  equations  (19), 
one  obtains,  after  some  manipulation, 
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(A  Ak+  B^)  - J - (AlAk-BlV 


/  ,  1  \ 

,  aR 

sin  (k  +  —  ^ 

)  — 

k4 


a-  "-j  / 

i  k  1  k 


k4 


sin(k-^)^-  c  .  3in  (k  +  R 

iAk +  B,  Bk>  — TT - +  (A^Ak  '  W  k  +  i 


(k4')f 


sin  (k  -t) 

Vk +  W  — rv 


sin(k+<0^ 


=  0, 
(20) 


(Vk  -  W  ~im  =  °’ 


(AA  +  W  kTT 


sin(k-t)“  sin(k+t)^ 

- r-V-  +  (A^Ak  -  B^)  - vTT -  =  ° 


k  +  l 


The  solutions  (19)  that  satisfy  (20)  and  have  real  coefficients  are  of 
the  form 

<°k  =  Ak  [Sin  k  (S  “  ^  )  ±  003  k  (  8  "  ~2~  )  1  ’ 

where  the  eigenvalues  are  given  by 

k  _  +  i_  ^  n  =  0,  ±1,  ±2,  ±3,  ... 


Whichever  sign  is  used  in  equation  (21),  it  must  be  used  for  all  of  the 

eigenfunctions  in  order  for  them  to  be  mutually  orthogonal.  At  the 

isolated  points  where  a  is  a  multiple  of  ir,  the  coefficients  Ak  and 

B  can  be  assigned  in  an  arbitrary  manner.  However,  since  we  are 
k 

concerned  here  with  .the  behavior  of  a  real  beam,  these  isolated  points 

are  of  no  particular  importance. 

Figure  2  shows  a  plot  of  the  three  lowest  eigenvalues  from 
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equation  (22)  vs.  beam  angle  a-  For  comparison  purposes  one  of  the 
classical  curves  from  equation  (6)  or  (7)  is  also  included.  The  ordinate 
of  this  graph  is  shown  positive  downwards  in  order  to  emphasize  that 
values  of  kR  less  than  1  correspond  to  a  positive  (closing)  moment 
and  values  of  kR  greater  than  1  to  a  negative  (opening)  moment.  ^  It 
can  be  seen  that  the  classical  curve  causes  the  moment  to  change  sign 
and  is  asymptotic  to  a  positive  value,  while  the  new  curves  are  asymp¬ 
totic  to  kR  =  1  which  corresponds  to  zero  buckling  moment.  The  small 
rectangles  on  this  figure  refer  to  experimental  results  which  will  be 
discussed  below. 

These  uniquely  determined  values  of  k  are  next  used  to  deter¬ 
mine  the  value  of  the  buckling  moment  M  from  equation  (4).  However, 
since  equation  (4)  is  quadratic  in  M,  extraneous  solutions  maybe  en¬ 
countered.  To  guard  against  these  extraneous  solutions,  we  note  that  the 
upper  curve  in  Figure  2  should, only  give  rise  to  positive  moments  and 
the  lower  curve  to  negative  moments.  In  addition,  the  magnitude  of  the 
buckling  moment  should  monotonically  decrease  with  increasing  a.  This 
arises  from  the  fact  that  a  curved  beam  can  buckle  into  the  true  buckled 
shape  of  a  shorter  curved  beam,  with  the  remaining  portion  having  only 
a  rigid  body  displacement.  This  "unbent"  portion  will  not  change  either 
the  potential  energy  or  the  boundary  conditions;  thus  the  longer  beam 

o.  \ 

This  behavior  would  have  been  reversed  had  k  =  -~  been  chosen 

as  the  fundamental  rigid  body  rotation. 
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could  buckle  at  the  same  buckling  moment  as  the  shorter  beam. 

Applying  these  considerations  to  the  buckling  moments  obtained 

from  Figure  2  and  equation  (4);  one  obtains  a  typical  buckling  moment 

curve  as  shown  in  Figure  3.  The  discontinuity  in  the  positive  buckling 

moment  curve  arises  from  the  fact  that  it  is  composed  of  portions  of 

two  different  curves.  The  lower  curve  passes  thru  zero  at  a  =  ir  and 

increases  until  it  reaches  a  maximum  at  a  ~  2tr.  The  upper  curve  would 

begin  increasing  for  a  greater  than  2n;  thus  the  decreasing  moment 

requirement  results  in  a  discontinuity  at  a  =  2ir,  where  the  buckling 

C  EIx 

moment  decreases  from  M  =  —  to  M  =  — — .  However,  this  discon- 

K.  K. 

tinuity  is  only  of  theoretical  interest,  since  any  real  curved  beam  must 
necessarily  have  a  less  than  2tt  by  some  small  amount.  In  passing, 
it  should  be  noted  that  these  buckling  moment  curves  are  the  minimum 
possible  for  lateral  buckling  of  a  curved  beam.  The  addition  of  any 
lateral  restraints  on  the  ends  of  the  beam  can  only  serve  to  increase  the 
buckling  moments. 

It  is  perhaps  helpful  to  plot  the  eigenfunctions  given  by  equation 
(21)  in  order  to  show  the  types  of  conditions  that  result  at  the  ends  of 
the  beam.  Figure  4  is  such  a  plot  for  the  case  of  buckling  under  a 
negative  moment.  Various  values  of  o:  are  represented  between  the 
limiting  cases  of  a  "  0  and  a  -  2rr.  Figures  5  and  6  present  the 
same  information  for  buckling  under  a  positive  moment  and  for  a  rigid 
body  rotation  about  a  diameter.  In  these  plots  it  is  helpful  to  remember 
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negat 


scrtteiaiiai 


that  k  represents  the  number  of  full  sine  waves  that  would  occur  if 
the  beam  angle  were  2tr.  These  curves  were  obtained  by  using  the  plus 
sign  in  equation  (21);  use  of  the  minus  sign  would  essentially  interchange 
the  ends.  Since  tp  can  be  related  to  v  by  equations  (1)  and  (2),  these 
curves  also  represent  the  deflection  of  the  beam  centerline  to  within  an 
arbitrary  rigid  body  motion.  Note  that,  as  would  be  expected,  the  eigen¬ 
functions  approach  the  same  shape  for  positive  or  negative  buckling  as 
the  beam  angle  approaches  zero.  Also,  the  lowest  classical  shape,  for 
either  positive  or  negative  bending  moments,  would  simply  be  half  of  a 
sine  wave  regardless  of  the  value  of  a. 
b.  Straight  beam 

One  can  obtain  the  buckling  moment  for  a  straight  beam  of  length 

L  by  taking  the  limit  of  the  curved  beam  equation  (22)  as  R  -♦  00  and 

2n7r 

aR  -  L.  This  gives  k  =  -jj-  which  is  twice  the  accepted  classical 
value  [4].  Because  of  this  disagreement,  it  is  necessary  .to  present 
the  analysis  although  it  is  similar  to  that  for  the  curved  beam. 

The  equations  for  the  straight  beam  can  be  readily  obtained 
from  equations  (1)  and  (2)  by  omitting  those  terms  which  contain  —  . 

In  this  manner  one  obtains 


i3 

d  (p  ,  .  a 

d#  +  k  V  = 


(23) 


where 


k3  = 


M 


El  C 
x 


(24) 
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19 


A.  A.  +  B.B 


sin  (k  -l)  — 


tk  Ik i  k-l 


sin  (k  -Yl)  - 

\W  "  Vk  )  k  +  l 


=  0  ,  (28c) 


siii  (k  -  t)  — 


sin  (k  +  £)  2 


(AA+  BA)  'T-T +  (AA -  ■  0  •  <28d> 


The  only  solution  to  the  system  (28)  which  has  zero  potential  energy 

for  the  trivial  solution  (0  ,  is 

o 

(0  =  B  , 

’o  o 

L  L 

( pk  =  Aj_  sin  k;{c.  -  -■)  +  cos  k(s--) 


(29) 


where 

k  =  ,  n  =  0,  ±1,  ±2,  .  .  .  .  (30) 

Li 

Note  that  since  the  A  and  B  are  arbitrary,  equations  (29)  and  (3  0) 

K  K 

contain  the  limit  of  equations  (21)  and  (22)  (as  R  -*  “  and  aR  -*  L) 
as  a  special  case. 


4.  EXPERIMENTAL  RESULTS 

In  order  to  verify  the  trend  of  the  results  given  by  equations 
(22)  and  (30),  a  short  experimental  program  was  undertaken.  Since 
there  is  a  relatively  wide  disagreement  between  these  results  and  the 
classical  values,  the  intent  of  the  experiment  was  not  to  provide  a 
rigorous  proof  of  the  new  results  but -rather  to  be  of  sufficient  accuracy 
to  show  the  disagreement  with  the  classical  results. 
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Before  the  buckling  moment  of  a  curved  beam  has  been  reached, 
the  stress  distribution  is  simply  that  of  the  well-known  plane  stress 
problem  of  the  pure  bending  of  a  curved  beam  [5].  Therefore  it  was 
decided  to  make  the  test  specimens  of  a  photoelastic  material  in  order 
to  have  visual  proof  that  the  loading  was  giving  the  correct  theoretical 
stress  distribution  befoxe  buckling  has  occurred.  The  specimens  were 
therefore  cut  from  photoelastic  material  PS-2B,  and  were  circular 
ring  sectors  of  5"  O.  D.  ,  2|r"  I.  D. ,  and  .081"  thick.  The  specimens 
were  then  epoxy  bonded  to  metal  strips  which  were  hinged  at  the  center 
of  the  circular  ring  sector.  Figure  7  shows  an  unstressed  full  ring  test 
specimen  in  the  transmission  polariscope  with  dark  field  illumination. 

In  this  case  the  metal  strip  was  actually  a  strap  hinge;  in  other  tests 
it  was  found  necessary  to  use  metal  strips  bolted  to  a  hinge,  so  that  the 
center  of  rotation  could,  be  adjusted  until  the  axisymmetric  plane  stress 
solution  was  obtained  upon  opening  or  closing  the  hinge.  Figure  8  shows 
the  same  specimen  after  pulling  upward  on  the  top  of  the  strap  hinge  to 
give  a  measured  opening  moment  of  -16.0  in- lb.  It  can  be  seen  that  the 
correct  axisymmetric  stress  distribution  has  been  obtained  except,  of 
course,  in  the  immediate  vicinity  of  the  bonded  ends.  The  first  dark 
fringe  is  just  at  the  inside  radius  of  the  specimen.  For  this  geometry 
and  material,  each  dark  fringe  corresponds  to  a  principle  stress  differ¬ 
ence  of  approximately  970  psi. 

Figure  9  shows  the  same  specimen  at  incipient  buckling,  and 
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Figure  7.  Stress  free  full  ring  test  specimen  in  the  transmission 
polariscope  with  dark  field  illumination. 


Figure  8.  Full  ring  specimen  with  an  opening  moment  of  -16.0 

in-lb.  The  stress  optical  coefficient  is  970  psi/fringe. 
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Figure  9-  Specimen  at  incipient  buckling  with  an  opening  moment 
of  approximately  -20.6  in-lb. 


Figure  10.  Post-buckled  specimen  with  an  opening  moment  of  -20.8 
in-lb.  The  dark  radial  bands  are  inflection  points. 
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Figure  10  shows  the  post-buckled  specimen  with  an  opening  moment  of 
-20.8  in-lb.  The  four  radial  dark  bands  inFigure  lOwere  experiment¬ 
ally  seen  to  correspond  to  the  inflection  points  where  tp  =  0.  This  post- 
buckled  shape  is,  therefore,  in  fairly  good  agreement  with  the  topmost 
theoretical  curve  shown  inFigure  4.  Note  thatthe  classical  post-buckled 
shape  of  equation  (6)  would  have  no  inflection  points  and  the  classical 
shape  of  equation  (7)  would  have  either  one  or  two  inflection  points,  de¬ 
pending  on  the  boundary  conditions.  The  strap  hinge  arrangement  was 
intended  to  be  loose  enough  laterally,  so  that  the  ends  would  behave  like 
the  "free"  edges  required  by  equation  (3).  However,  it  is  possible  that 
St.  Venant's  principle  applies  to  this  problem,  since  it  was  observed 
that  a  small  deflection  of  the  ends  did  not  seem  to  contribute  to  the  lateral 
buckling,  in  contrast  to  the  problem  of  buckling  under  an  axial  force. 

Figures  11  and  12  show  the  pre-  and  post-buckled  stress  dis¬ 
tributions  for  a  three-quarter  ring  specimen  which  buckled  under  an 
opening  moment  of  -27.7  in-lb.  The  three  inflection  points  in  Figure  12 
can  be  seen  to  be  in  reasonable  agreement  with  the  second  theoretical 
curve  of  Figure  4.  Figures  13  and  14  show  the  same  information  for  a 
half  ring  specimen  which  buckled  at  -41.8  in- lb.  Some  difficulty  was 
encountered  in  maintaining  the  plane  stress  solution  for  this  specimen, 
and  any  slight  departure  from  the  plane  stress  solution  caused  a  con¬ 
siderable  reduction  in  the  buckling  moment.  In  addition,  the  lateral 
restraints  had  to  be  maintained  so  loosely  that  the  buckling  process 


Figure  13.  Half  ring  specimen  with  an  opening  moment 
of  -21.9  in-lb.. 


Figure  14.  Post-buckled  specimen  which  buckled  at  an 
opening  moment  of  -41.8  in- lb. 


itself  caused  a  considerable  change  in  the  stress  distribution.  Con¬ 
sequently,  the  buckled  shape  in  Figure  14  is  only  in  rough  agreement 
with  the  third  theoretical  curve  of  Figure  4. 

Data  was  also  obtained  for  buckling  under  a  closing  moment. 
However,  because  of  a  shortage  of  photoelastic  material,  the  larger 
specimens  were  simply  cut  down  to  make  the  shorter  specimens;  and, 
due  to  an  unfortunate  oversight,  no  data  was  obtained  on  the  approxi¬ 
mately  full  ring  specimen.  The  three-quarter  ring  specimen,  however, 
is  shown  in  Figures  15  and  16  for  pre-  and  post-buckling  under  a  closing 
moment.  The  absence  of  an  inflection  point  is  in  agreement  with  the 
second  theoretical  curve  of  Figure  5.  Figure  17  shows  the  half  ring 
specimen  with  a  closing  moment  of  approximately  49  in-lb.  A  slight 
increase  in  this  moment  appeared  to  cause  the  specimen  to  simultan¬ 
eously  buckle  and  fracture.  Because  of  this  tensile  strength  limitation, 
no  data  on  shorter  specimens  was  obtained. 

The  experimentally  obtained  buckling  moments  were  used  to 
calculate  values  of  k  for  these  five  specimens.  These  values  are 
indicated  by  the  small  squares  on  Figure  2.  In  addition,  the  buckling 
moments  themselves  are  shown  on  Figure  3.  The  amount  of  agreement 
between  the  data  points  and  the  theoretical  curves  should  be  considered 
fortuituous,  since  there  were  uncertainties  in  both  the  experimental 
procedure  and  the  material  properties  of  the  specimens. 

Figures  18  and  19  show  the  pre-  and  post-buckled  stress 
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Figure  17.  Half  ring  specimen  immediately  before  failure  with 
a  closing  moment  of  +49.0  in- lb. 
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Figure  ,8.  Straight  beam  specimen  with  a  calculated 
applied  moment  of  20.0  in-lb. 


Figure  1 9.  Post-buckled 


specimen. 
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distributions  for  an  initially  straight  beam  of  the  same  material.  The 
beam  dimensions  were  l\u  long  by  0.8"  high  and  .081"  thick.  The 
specimen  was  epoxy  bonded  to  long  transverse  metal  strips;  the  right 
hand  strip  was  clamped  to  the  testing  apparatus  and  the  left  hand  strip 
remained  free.  The  bending  moment  was  then  applied  by  hand  with  care 
being  taken  to  obtain  the  plane  stress  solution  before  buckling  had  occured. 
The  applied  bending  moment  was,  therefore,  calculated  from  the  observed 
stress  pattern,  using  the  dimensions  of  the  beam  and  the  stress  optical 
coefficient  of  the  material. 

The  specimen  remained  straight  above  the  classical  buckling 
moment  and  had  a  tendency  to  return  to  the  straight  form  when  a  very 
small  buckle  was  manually  introduced.  The  specimen  finally  buckled 
when  the  calculated  applied  moment  reached  a  value  of  20  in-lb,  which 
is  almost  exactly  twice  the  lowest  classical  value  for  these  dimensions; 
this  is  in  agreement  with  equation  (30).  Note  the  similarity  of  the 
buckled  shape  of  Figure  19  to  the  bottommost  theoretical  curve  of  either 
Figure  4  or  Figure  This  is  possibly  due  to  the  fact  that  ^he  straight 
beam  is  actually  significantly  curved  before  buckling  finally  occurs. 

5,  SUMMARY 

The  basic  differential  equation  for  the  lateral  buckling  of  a  beam 
subjected  to  pure  bending  is  derived  under  the  assumption  that  there  are 
no  lateral  forces  acting  on  the  ends  of  the  beam.  Therefore,  the 
classical  approach  of  assigning  relatively  arbitrary  boundary  conditions 
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will,  in  general,  violate  this  assumption,  and  there  will  be  shearing 
forces  on  the  boundaries  which  were  not  included  in  the  differential 
equation.  In  contrast,  the  method  presented  in  this  paper  insures  that 
the  principle  of  stationary  potential  energy  (i.  e.  ,  virtual  work)  is  satis¬ 
fied  and  the  resulting  orthogonality  conditions  lead  to  a  unique  set  of 
boundary  conditions.  The  experimental  data  agrees  with  these  new 
results  rather  than  the  classical  results  —  in  magnitude  of  the  buckling 
moment,  as  well  as  in  the  post-buckled  shape  of  the  beam. 


APPENDIX 


1 .  Orthogonality  Relationships  for  a  Straight  Beam  Subjected  to  an 
Axial  Force. 

For  an  initially  straight  beam  subjected  to  an  axial  force  P 
whose  direction  of  action  remains  fixed  in  space,  the  equilibrium 
equation  for  the  deflection  v  of  the  beam  is 


a) 


Since  the  force  P  is  the  only  agent  doing  work  oh  the  beam,  the 
potential  energy  expression  for  the  beam,  in  the  case  of  an  inextensible 
centerline,  is  given  by 


V  = 


P  (WdvA3 
ds  J 


ds 


(2) 


We  now  consider  that  we  have  a  complete  set  of  solutions  v  to 

K 

equation  (1)  which  satisfy  the  same  natural  boundary  conditions 
giving  rise  to  (I)  and  (2);  i.  e.  ,  either  the  moment  or  the  slope 
is  zero  on  the  boundaries,  and,  either  the  shear  or  the  displacement 
is  zero  on  the  boundaries.  We  now  apply  the  principle  of  stationary 
potential  energy  at  the  eigenstates  v^  and  v^.  Taking  a  varied  state 
about  the  true  state  v^  in  the  form  v*  =  v^  +  e  v^,  we  obtain  from 
equation  (2), 


a  v  (vi,;) 
a  e 


o  = 


ei  r 


L  d"  v,  d  v 


P,  L  dv,  dv. 


V  .  A  .  a-#  V*  V  ,  V  » 

)(t>-  »> 
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Taking  a  varied  state  about  the  true  state  v  in  the  form  v**  =  v.  +  e  v,  > 

v  1/  K 

we  obtain 


bV  (y*'c) 

be 


L  d^v,  d2v. 


P,  L  dv,  dv, 
l  P  /  k\  / 


*  <*  ■  ?  J0  -  tJ  (ir)(ir>-  «> 

€  "*  1/  U  U 


Therefore,  if  E  7^  P,  ,  we  obtain  from  equations  (3)  and  (4)  the 

following  orthogonality  relations  for  the  eigenfunctions: 

L  d2  v,  d2  v 


L  dv,  dv. 

I  (<rk)  (it4) 


ds  = 


(5) 
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This  paper  presents  a  new  analysis  of  the  classical  problem  of  the 
lateral  buckling  under  pure  bending  of  a  curved  beam  with  circular  axis. 
Although  a  solution  to  this  problem,  and  to  the  corresponding  straight  beam 
problem  is  well-accepted  in  the  literature,  the  present  analysis  gives  buck¬ 
ling  moments  which  are  significantly  different  than  the  usually  accepted  results. 
The  difference  arises  from  an  application  of  new  boundary  conditions  that  are 
shown  io  be  the  only  boundary  conditions  consistent  with  the  equilibrium  of  the 
buckled  beam  under  a  pure  bending  moment.  This  new  analysis  provides 
generally  higher  bending  moments.  A  short  experimental  program  was  per¬ 
formed,  and  the  results  tend  to  confirm  the  new  analysis. 
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